The synergy between the PscC subunits for electron transfer to the P840 special pair in Chlorobaculum tepidum.

The primary photochemical reaction of photosynthesis in green sulfur bacteria occurs in the homodimer PscA core proteins by a special chlorophyll pair. The light induced excited state of the special pair producing P840+ is rapidly reduced by electron transfer from one of the two PscC subunits. Molecular dynamics (MD) simulations are combined with bioinformatic tools herein to provide structural and dynamic insight into the complex between the two PscA core proteins and the two PscC subunits. The microscopic dynamic model involves extensive sampling at atomic resolution and at a cumulative time-scale of 22µs and reveals well defined protein-protein interactions. The membrane complex is composed of the two PscA and the two PscC subunits and macroscopic connections are revealed within a putative electron transfer pathway from the PscC subunit to the special pair P840 located within the PscA subunits. Our results provide a structural basis for understanding the electron transport to the homodimer RC of the green sulfur bacteria. The MD based approach can provide the basis to further probe the PscA-PscC complex dynamics and observe electron transfer therein at the quantum level. Furthermore, the transmembrane helices of the different PscC subunits exert distinct dynamics in the complex.

Electric Field Susceptibility of Chlorophyll c Leads to Unexpected Excitation Dynamics in the Major Light-Harvesting Complex of Diatoms.

Diatoms are one of the most abundant photosynthetic organisms on earth and contribute largely to atmospheric oxygen production. They contain fucoxanthin and chlorophyll-a/c binding proteins (FCPs) as light-harvesting complexes with a remarkable adaptation to the fluctuating light on ocean surfaces. To understand the basis of the photosynthetic process in diatoms, the excitation energy funneling within FCPs must be probed. A state-of-the-art multiscale analysis within a quantum mechanics/molecular mechanics framework has been employed. To this end, the chlorophyll (Chl) excitation energies within the FCP complex from the diatom Phaeodactylum tricornutum have been determined. The Chl-c excitation energies were found to be 5-fold more susceptible to electric fields than those of Chl-a pigments and thus are significantly lower in FCP than in organic solvents. This finding challenges the general belief that the excitation energy of Chl-c is always higher than that of Chl-a in FCP proteins and reveals that Chl-c molecules are much more sensitive to electric fields within protein scaffolds than in Chl-a pigments. The analysis of the linear absorption spectrum and the two-dimensional electronic spectra of the FCP complex strongly supports these findings and allows us to study the excitation transfer within the FCP complex.

Engineering salt-tolerant Cas12f1 variants for gene-editing applications.

CRISPR has revolutionized the field of genome editing in life sciences by serving as a versatile and state-of-the-art tool. Cas12f1 is a small nuclease of the bacterial immunity CRISPR system with an ideal size for cellular delivery, in contrast to CRISPR-associated (Cas) proteins like Cas9 or Cas12. However, Cas12f1 works best at low salt concentrations. In this study, we find that the plasticity of certain Cas12f1 regions (K196-Y202 and I452-L515) is negatively affected by increased salt concentrations. On this line, key protein domains (REC1, WED, Nuc, lid) that are involved in the DNA-target recognition and the activation of the catalytic RuvC domain are in turn also affected. We suggest that salt concentration should be taken in to consideration for activity assessments of Cas engineered variants, especially if the mutations are on the protospacer adjacent motif interacting domain. The results can be exploited for the engineering of Cas variants and the assessment of their activity at varying salt concentrations. We propose that the K198Q mutation can restore at great degree the compromised plasticity and could potentially lead to salt-tolerant Cas12f1 variants. The methodology can be also employed for the study of biomolecules in terms of their salinity tolerance.Communicated by Ramaswamy H. Sarma.

Corrigendum to "Recognition motifs for importin 4 [(L)PPRS(G/P)P] and importin 5 [KP(K/Y)LV] binding, identified by bio-informatic simulation and experimental in vitro validation" [Comput Struct Biotechnol J 20 (2022) 5952-5961].

[This corrects the article DOI: 10.1016/j.csbj.2022.10.015.].

A perspective on the major light-harvesting complex dynamics under the effect of pH, salts, and the photoprotective PsbS protein.

The photosynthetic apparatus is a highly modular assembly of large pigment-binding proteins. Complexes called antennae can capture the sunlight and direct it from the periphery of two Photosystems (I, II) to the core reaction centers, where it is converted into chemical energy. The apparatus must cope with the natural light fluctuations that can become detrimental to the viability of the photosynthetic organism. Here we present an atomic scale view of the photoprotective mechanism that is activated on this line of defense by several photosynthetic organisms to avoid overexcitation upon excess illumination. We provide a complete macroscopic to microscopic picture with specific details on the conformations of the major antenna of Photosystem II that could be associated with the switch from the light-harvesting to the photoprotective state. This is achieved by combining insight from both experiments and all-atom simulations from our group and the literature in a perspective article.

Recognition motifs for importin 4 [(L)PPRS(G/P)P] and importin 5 [KP(K/Y)LV] binding, identified by bio-informatic simulation and experimental in vitro validation.

Nuclear translocation of large proteins is mediated through karyopherins, carrier proteins recognizing specific motifs of cargo proteins, known as nuclear localization signals (NLS). However, only few NLS signals have been reported until now. In the present work, NLS signals for Importins 4 and 5 were identified through an unsupervised in silico approach, followed by experimental in vitro validation. The sequences LPPRS(G/P)P and KP(K/Y)LV were identified and are proposed as recognition motifs for Importins 4 and 5 binding, respectively. They are involved in the trafficking of important proteins into the nucleus. These sequences were validated in the breast cancer cell line T47D, which expresses both Importins 4 and 5. Elucidating the complex relationships of the nuclear transporters and their cargo proteins is very important in better understanding the mechanism of nuclear transport of proteins and laying the foundation for the development of novel therapeutics, targeting specific importins.

From Traditional Ethnopharmacology to Modern Natural Drug Discovery: A Methodology Discussion and Specific Examples.

Ethnopharmacology, through the description of the beneficial effects of plants, has provided an early framework for the therapeutic use of natural compounds. Natural products, either in their native form or after crude extraction of their active ingredients, have long been used by different populations and explored as invaluable sources for drug design. The transition from traditional ethnopharmacology to drug discovery has followed a straightforward path, assisted by the evolution of isolation and characterization methods, the increase in computational power, and the development of specific chemoinformatic methods. The deriving extensive exploitation of the natural product chemical space has led to the discovery of novel compounds with pharmaceutical properties, although this was not followed by an analogous increase in novel drugs. In this work, we discuss the evolution of ideas and methods, from traditional ethnopharmacology to in silico drug discovery, applied to natural products. We point out that, in the past, the starting point was the plant itself, identified by sustained ethnopharmacological research, with the active compound deriving after extensive analysis and testing. In contrast, in recent years, the active substance has been pinpointed by computational methods (in silico docking and molecular dynamics, network pharmacology), followed by the identification of the plant(s) containing the active ingredient, identified by existing or putative ethnopharmacological information. We further stress the potential pitfalls of recent in silico methods and discuss the absolute need for in vitro and in vivo validation as an absolute requirement. Finally, we present our contribution to natural products' drug discovery by discussing specific examples, applying the whole continuum of this rapidly evolving field. In detail, we report the isolation of novel antiviral compounds, based on natural products active against influenza and SARS-CoV-2 and novel substances active on a specific GPCR, OXER1.

Deciphering the QR Code of the CRISPR-Cas9 System: Synergy between Gln768 (Q) and Arg976 (R).

Markov state models (MSMs) and machine learning (ML) algorithms can extrapolate the long-time-scale behavior of large biomolecules from molecular dynamics (MD) trajectories. In this study, an MD-MSM-ML scheme has been applied to probe the large endonuclease (Cas9) in the bacterial adaptive immunity CRISPR-Cas9 system. CRISPR has become a programmable and state-of-the-art powerful genome editing tool that has already revolutionized life sciences. CRISPR-Cas9 is programmed to process specific DNA sequences in the genome. However, human/biomedical applications are compromised by off-target DNA damage. Characterization of Cas9 at the structural and biophysical levels is a prerequisite for the development of efficient and high-fidelity Cas9 variants. The Cas9 wild type and two variants (R63A-R66A-R70A, R69A-R71A-R74A-R78A) are studied herein. The configurational space of Cas9 is provided with a focus on the conformations of the side chains of two residues (Gln768 and Arg976). A model for the synergy between those two residues is proposed. The results are discussed within the context of experimental literature. The results and methodology can be exploited for the study of large biomolecules in general and for the engineering of more efficient and safer Cas9 variants for applications.

Natural Polyphenols Inhibit the Dimerization of the SARS-CoV-2 Main Protease: The Case of Fortunellin and Its Structural Analogs.

3CL-Pro is the SARS-CoV-2 main protease (MPro). It acts as a homodimer to cleave the large polyprotein 1ab transcript into proteins that are necessary for viral growth and replication. 3CL-Pro has been one of the most studied SARS-CoV-2 proteins and a main target of therapeutics. A number of drug candidates have been reported, including natural products. Here, we employ elaborate computational methods to explore the dimerization of the 3CL-Pro protein, and we formulate a computational context to identify potential inhibitors of this process. We report that fortunellin (acacetin 7-O-neohesperidoside), a natural flavonoid O-glycoside, and its structural analogs are potent inhibitors of 3CL-Pro dimerization, inhibiting viral plaque formation in vitro. We thus propose a novel basis for the search of pharmaceuticals as well as dietary supplements in the fight against SARS-CoV-2 and COVID-19.

Robust Strategy for Photoprotection in the Light-Harvesting Antenna of Diatoms: A Molecular Dynamics Study.

Diatoms generate a large portion of the oxygen produced on earth due to their exceptional light-harvesting properties involving fucoxanthin and chlorophyll-binding proteins (FCP). At the same time, an efficient adaptation of these complexes to fluctuating light conditions is necessary to protect the diatoms against photodamage. So far, structural and dynamic data for the interaction between FCP and the photoprotective LHCX family of proteins in diatoms are lacking. In this computational study, we provide a structural basis for a remarkable pH-dependent adaptation at the molecular level. Upon binding of the LHCX1 protein to the FCP complex together with a change in pH, conformational changes within the FCP protein result in a variation of the electronic coupling in a specific chlorophyll-fucoxanthin pair, leading to a change in the exciton transfer rate by almost an order of magnitude. A common strategy for photoprotection between diatoms and higher plants is identified and discussed.

Time-dependent atomistic simulations of the CP29 light-harvesting complex.

Light harvesting as the first step in photosynthesis is of prime importance for life on earth. For a theoretical description of photochemical processes during light harvesting, spectral densities are key quantities. They serve as input functions for modeling the excitation energy transfer dynamics and spectroscopic properties. Herein, a recently developed procedure is applied to determine the spectral densities of the pigments in the minor antenna complex CP29 of photosystem II, which has recently gained attention because of its active role in non-photochemical quenching processes in higher plants. To this end, the density functional-based tight binding (DFTB) method has been employed to enable simulation of the ground state dynamics in a quantum-mechanics/molecular mechanics (QM/MM) scheme for each chlorophyll pigment. Subsequently, the time-dependent extension of the long-range corrected DFTB approach has been used to obtain the excitation energy fluctuations along the ground-state trajectories also in a QM/MM setting. From these results, the spectral densities have been determined and compared for different force fields and to spectral densities from other light-harvesting complexes. In addition, time-dependent and time-independent excitonic Hamiltonians of the system have been constructed and applied to the determination of absorption spectra as well as exciton dynamics.

p-cymene impairs SARS-CoV-2 and Influenza A (H1N1) viral replication: In silico predicted interaction with SARS-CoV-2 nucleocapsid protein and H1N1 nucleoprotein.

Therapeutic regimens for the COVID-19 pandemics remain unmet. In this line, repurposing of existing drugs against known or predicted SARS-CoV-2 protein actions have been advanced, while natural products have also been tested. Here, we propose that p-cymene, a natural monoterpene, can act as a potential novel agent for the treatment of SARS-CoV-2-induced COVID-19 and other RNA-virus-induced diseases (influenza, rabies, Ebola). We show by extensive molecular simulations that SARS-CoV-2 C-terminal structured domain contains a nuclear localization signal (NLS), like SARS-CoV, on which p-cymene binds with low micromolar affinity, impairing nuclear translocation of this protein and inhibiting viral replication, as verified by preliminary in vitro experiments. A similar mechanism may occur in other RNA-viruses (influenza, rabies and Ebola), also verified in vitro for influenza, by interaction of p-cymene with viral nucleoproteins, and structural modification of their NLS site, weakening its interaction with importin A. This common mechanism of action renders therefore p-cymene as a possible antiviral, alone, or in combination with other agents, in a broad spectrum of RNA viruses, from SARS-CoV-2 to influenza A infections.

Hydrogen gas as a central on-off functional switch of reversible metabolic arrest - New perspectives for biotechnological applications.

Metabolism is the sum of all chemical reactions that sustain life. There is an ongoing effort to control metabolic rate, which correlates with the maximum lifespan potential and constitutes one of the oldest scientific questions. Herein, we report on the complete reversible arrest of cellular metabolism and cell growth in a series of organisms, from microalgae to yeast upon exposure to a 100 % hydrogen atmosphere. We also report a tolerance of the microalgae under these conditions against extreme stress conditions, like high salt concentrations. The addition of oxygen or air almost completely restores the metabolic rate and cell growth. Molecular dynamics simulations are employed to decipher this phenomenon at atomic scale. Various proteins, including photosynthetic and respiratory complexes (LHCII, cytochrome c5) are probed in the interaction with hydrogen. Exposure to hydrogen, as opposed to oxygen, decreases the fluctuations of protein residues indicating thermostability. According to the above mechanism, an absolute hydrogen atmosphere can preserve biological products (e.g. fruits) for a long time without consuming any energy. By combining biological, chemical and computational methods, in this research we provide the basis for future innovative studies and advances in the field of biotechnology.

Multiscale QM/MM molecular dynamics simulations of the trimeric major light-harvesting complex II.

Photosynthetic processes are driven by sunlight. Too little of it and the photosynthetic machinery cannot produce the reductive power to drive the anabolic pathways. Too much sunlight and the machinery can get damaged. In higher plants, the major Light-Harvesting Complex (LHCII) efficiently absorbs the light energy, but can also dissipate it when in excess (quenching). In order to study the dynamics related to the quenching process but also the exciton dynamics in general, one needs to accurately determine the so-called spectral density which describes the coupling between the relevant pigment modes and the environmental degrees of freedom. To this end, Born-Oppenheimer molecular dynamics simulations in a quantum mechanics/molecular mechanics (QM/MM) fashion utilizing the density functional based tight binding (DFTB) method have been performed for the ground state dynamics. Subsequently, the time-dependent extension of the long-range-corrected DFTB scheme has been employed for the excited state calculations of the individual chlorophyll-a molecules in the LHCII complex. The analysis of this data resulted in spectral densities showing an astonishing agreement with the experimental counterpart in this rather large system. This consistency with an experimental observable also supports the accuracy, robustness, and reliability of the present multi-scale scheme. To the best of our knowledge, this is the first theoretical attempt on this large complex system is ever made to accurately simulate the spectral density. In addition, the resulting spectral densities and site energies were used to determine the exciton transfer rate within a special pigment pair consisting of a chlorophyll-a and a carotenoid molecule which is assumed to play a role in the balance between the light harvesting and quenching modes.

Trivial Excitation Energy Transfer to Carotenoids Is an Unlikely Mechanism for Non-photochemical Quenching in LHCII.

Higher plants defend themselves from bursts of intense light via the mechanism of Non-Photochemical Quenching (NPQ). It involves the Photosystem II (PSII) antenna protein (LHCII) adopting a conformation that favors excitation quenching. In recent years several structural models have suggested that quenching proceeds via energy transfer to the optically forbidden and short-lived S 1 states of a carotenoid. It was proposed that this pathway was controlled by subtle changes in the relative orientation of a small number of pigments. However, quantum chemical calculations of S 1 properties are not trivial and therefore its energy, oscillator strength and lifetime are treated as rather loose parameters. Moreover, the models were based either on a single LHCII crystal structure or Molecular Dynamics (MD) trajectories about a single minimum. Here we try and address these limitations by parameterizing the vibronic structure and relaxation dynamics of lutein in terms of observable quantities, namely its linear absorption (LA), transient absorption (TA) and two-photon excitation (TPE) spectra. We also analyze a number of minima taken from an exhaustive meta-dynamical search of the LHCII free energy surface. We show that trivial, Coulomb-mediated energy transfer to S 1 is an unlikely quenching mechanism, with pigment movements insufficiently pronounced to switch the system between quenched and unquenched states. Modulation of S 1 energy level as a quenching switch is similarly unlikely. Moreover, the quenching predicted by previous models is possibly an artifact of quantum chemical over-estimation of S 1 oscillator strength and the real mechanism likely involves short-range interaction and/or non-trivial inter-molecular states.

The conformational phase space of the photoprotective switch in the major light harvesting complex II.

Transitions between protein states are triggered by external stimuli. This knowledge leads to the control of protein function. Herein, we report a large scale (90 µs) study on the conformational space of the major light harvesting complex II, based on a comprehensive array of external stimuli.

Structure and Dynamics of a Thermostable Alcohol Dehydrogenase from the Antarctic Psychrophile Moraxella sp. TAE123.

The structure of a recombinant (His-tagged at C-terminus) alcohol dehydrogenase (MoADH) from the cold-adapted bacterium Moraxella sp. TAE123 has been refined with X-ray diffraction data extending to 1.9 Å resolution. The enzyme assumes a homo-tetrameric structure. Each subunit comprises two distinct structural domains: the catalytic domain (residues 1-150 and 288-340/345) and the nucleotide-binding domain (residues 151-287). There are two Zn2+ ions in each protein subunit. Two additional zinc ions have been found in the crystal structure between symmetry-related subunits. The structure has been compared with those of homologous enzymes from Geobacillus stearothermophilus (GsADH), Escherichia coli (EcADH), and Thermus sp. ATN1 (ThADH) that thrive in environments of diverse temperatures. Unexpectedly, MoADH has been found active from 10 to at least 53 °C and unfolds at 89 °C according to circular dichroism spectropolarimetry data. MoADH with substrate ethanol exhibits a small value of activation enthalpy ΔH ‡ of 30 kJ mol-1. Molecular dynamics simulations for single subunits of the closely homologous enzymes MoADH and GsADH performed at 280, 310, and 340 K showed enhanced wide-ranging mobility of MoADH at high temperatures and generally lower but more distinct and localized mobility for GsADH. Principal component analysis of the fluctuations of both ADHs resulted in a prominent open-close transition of the structural domains mainly at 280 K for MoADH and 340 K for GsADH. In conclusion, MoADH is a very thermostable, cold-adapted enzyme and the small value of activation enthalpy allows the enzyme to function adequately at low temperatures.

Structural Basis for Allosteric Regulation in the Major Antenna Trimer of Photosystem II.

The allosteric regulation of protein function proves important in many life-sustaining processes. In plant photosynthesis, LHCII, the major antenna complex of Photosystem II, employs a delicate switch between light harvesting and photoprotective modes. The switch is triggered by an enlarged pH gradient (ΔpH) across the thylakoid membranes. Using molecular simulations and quantum calculations, we show that ΔpH can tune the light-harvesting potential of the antenna via allosteric regulation of the excitonic coupling in chlorophyll-carotenoid pairs. To this end, we propose how the LHCII excited state lifetime is coupled to the environmental conditions. In line with experimental findings, our theoretical model provides crucial evidence toward the elucidation of the photoprotective switch of higher plants at an all-atom resolution.

Fine tuning of the photosystem II major antenna mobility within the thylakoid membrane of higher plants.

Depending on the amount of light, the photosystem II (PSII) antennae or Light Harvesting Complexes (LHCII) switch between two states within the thylakoid membranes of higher plants, i.e., a light-harvesting and a photoprotective mode. This switch is co-regulated by a pH gradient (ΔpH) across the membrane and the interaction with the PSII subunit S (PsbS) that is proposed to induce LHCII aggregation. Herein we employ all-atom and coarse-grained molecular simulations of the major LHCII trimer at low and excess ΔpH, as well as in complexation with PsbS within a native thylakoid membrane model. Our results demonstrate the aggregation potential of LHCII and, consistent with the experimental literature, reveal the role of PsbS at atomic resolution. PsbS alters the LHCII-thylakoid lipid interactions and restores the LHCII mobility that is lost in the transition to photoprotective conditions (low lumenal pH). In agreement with this finding, diffusion of the integral membrane protein LHCII is dependent on both, electrostatic interactions and hydrophobic mismatch, while it does not obey the Saffman-Delbrück diffusion model.

Protein-protein interactions within photosystem II under photoprotection: the synergy between CP29 minor antenna, subunit S (PsbS) and zeaxanthin at all-atom resolution.

The assembly and disassembly of protein complexes within cells are crucial life-sustaining processes. In photosystem II (PSII) of higher plants, there is a delicate yet obscure balance between light harvesting and photo-protection under fluctuating light conditions, that involves protein-protein complexes. Recent breakthroughs in molecular dynamics (MD) simulations are combined with new approaches herein to provide structural and energetic insight into such a complex between the CP29 minor antenna and the PSII subunit S (PsbS). The microscopic model involves extensive sampling of bound and dissociated states at atomic resolution in the presence of photo-protective zeaxanthin (Zea), and reveals well defined protein-protein cross-sections. The complex is placed within PSII, and macroscopic connections are emerging (PsbS-CP29-CP24-CP47) along the energy transfer pathways from the antenna to the PSII core. These connections explain macroscopic observations in the literature, while the previously obscured atomic scale details are now revealed. The implications of these findings are discussed in the context of the Non-Photochemical Quenching (NPQ) of chlorophyll fluorescence, the down-regulatory mechanism of photosynthesis, that enables the protection of PSII against excess excitation load. Zea is found at the PsbS-CP29 cross-section and a pH-dependent equilibrium between PsbS dimer/monomers and the PsbS-CP29 dissociation/association is identified as the target for engineering tolerant plants with increased crop and biomass yields. Finally, the new MD based approaches can be used to probe protein-protein interactions in general, and the PSII structure provided can initiate large scale molecular simulations of the photosynthetic apparatus, under NPQ conditions.